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بهبود روش اجزای محدود سط یافته در آنالیز مسائل ترک با تولید مجدد شبکه * 
کاسرآن مزب (6 تک رتاش ۳ 


چکیه۵ روش اجزای محدود بسطیافته یکی از تکنیک‌های بسیار موثر برای مدل‌سازی مسائل دارای ترک می‌باشد. در این روش باتوجه به 
استفاده از توابع غنی‌سازی برای مد لکردن ناپیوستگی میدان جابه‌جایی» دیگر ترک به‌صورت یک ماهیت هندسی مدل نمی‌شود و نیازی نمی‌باشد 
که شیک اجزای محدود در هر مرحله مسیر ترک را دنبال کند. ولی مشکل اساسی در این روش عدم کنترل برروی میزان عطای حل عددی 
می‌باشد. به این منظور به‌کمک روش تولید مجدد شبکه نطابقی, در مواردی که حطای حل بیش از حد قابل قبول است, شبکهُ اجزای محدود 
مجددا تولید می‌شود و مسئله با شبکه جدید حل می‌شود تا بتوان به دقت مطلوب دست یافت. این تولید مجدد شبکه کاملا مستقل از دنبال 
کردن سیر ترک می‌باشد و ممکن است در چند مرحله از رشد ترک, مسئله با یک شبکه یکسان حل شود. الگوریتم ذکرشده باعث می‌شود که 
بدون مشکلات دنبال کردن مسیر ترک, دقت حل نیز در حد مطلوبی باقی بماند. پارامتر ضریب شدت تنش که نق شسکلیدی در نعیین مسیر رشد 


ترک دارد, با کمک روش ت رکیبی محاسبه شده‌است و بهبود قابل توجهی در دقت محاسبهٌ عددی این پارامتر مشاهده شده‌است. 


واژه‌های کلیدی اجزای محدود بسطیافته آنالیز مسائل ترک» اجزای محدود تطابقی. توابع غنی‌سازی, اصلاح مش. 
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*تاریخ دریافت مقاله ۹۵/۹/۵ تاریخ پذیرش آن ٩۱/4/۱۷‏ می‌باشد. 
(۱) کارشناس ارشد سازه, مهندسی عمران, دانشگاه شاهد. تهران. 
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مقد مه 


حل مسائل ترک به‌دلیل وجود ناپیوستگی در میدان 
جابه‌جایی و تکین بودن میدان تتش در نوک ترک همواره 
جزو مسائل چالش‌برانگیز در رشته‌های مختلف مهندسی 
بوده‌است. به‌حصوص در استفاده از روش‌های عددی. 
روش‌هایی که در مسائل بدون ترک به‌خوبی جواب مسئله 
را مدل می‌کنند. در مسائل دارای ترک جواب‌های واگرا 
و کاملا دور از جواب تحلیلی مسئله را می‌دهند. بنابراین 
برای بهبود روش‌های عددی در تحلیل مسائل ترک 
تکنیک‌های مختلفی توسط محققان مختلف پيشنهاد 
گشته‌است. یکی از این تکنیک‌ها روش اجزای محدود 
بسطيافته می‌باشد که ایدة اصلی این روش برمبنای اصل 
تقسیم‌بندی جزء واحد است که توسط بابوشکا و ملنک 
ارائه گردید [1], 

به‌موازات این روش. تکنیک اجزای محدود 
تطابقی یا همان روش اصلاح مش برای مدل‌سازی دقیق 
رفتار مسائل دارای ترک مورداستفاده قرار گرفته‌است. این 
روش نیز سعی دارد با ارائٌ اصلاح شبکه. خطاهای ایجاد 
فا که نت یی ید فا تتارش اش بات کاعی 
دهد. 

روش اجزای محدود بسطیافته یکی از تکنیک‌های 
موفق درزمينة آنالیز مسائل ترک می‌باشد. مهم‌ترین مزیتی 
که این روش نسبت به روش اجزای محدود معمولی دارد 
این است که دیگر نیازی به مش‌بندی مجدد در مسئله به 
هنگام رشد ترک نیست. در مسائلی که بااستفاده از روش 
اجزای محدود معمولی آنالیز می‌شوند در هر مرحله از 
رشد ترک هندسه مسئله کاملا تغییر می‌پابد چرا که 
بایستی خود ترک در این مسائل مستقیما مدل گردد و 
همین موضوع باعث برهم خوردن هندسة مسئله درزمانی 
هی اتیواق کر کوش من کی [2: 

اما در روش اجزای محدود بسطیافته با تغییر توابع 
شکل و افزودن جملاتی به معادلة جابه‌جایی. ترک یا 
ناپیوستگی مدل می‌شود و نیازی به مش‌بندی مجدد در 
هر مرحله از رشد ترک وجود ندارد. شکل () تفاوت 
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بهبود روش اجزای محدود بسط‌یافته د رآنالیز. 


مش‌بندی در دو روش اجزای محدود معمولی و بسط یافته 
را به تصویر می‌کشد [3]. 

برای مدل‌سازی ناپیوستگی در میدان جابه‌جایی ترک از 
افزودن جملاتی به تابع جابه‌جایی در روش اجزای 
محدود کلاسیک استفاده می‌شود. این جملات توابع 
غنی‌سازی نام دارند. توابع غنی‌سازی مورداستفاده در 
روش اجزای محدود بسطیافته به دو دسته تقسیم 
می‌شوند. دستهٌ اول توابع برای مدل‌سازی المان‌هایی 
می‌باشد که ترک از آنها عبور کرده‌است و دسته دوم برای 
المان‌هایی می‌باشد که نوک ترک در آنها واقع شده‌است. 


۷« 
شکل ۱ تطبیق و عدم تطبیق مش‌بندی در دو روش اجزای 
محدود کلاسیک و بسط یافته 


بلیچکو [4] برای ترک الاستیک و برای مسائل ایزوتروپ 


رابطهٌ (۱) را به‌عنوان رابطهٌ نهایی برای غنی‌سازی معرفی 


کرده‌است. 
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در روابط فوق لا مجموعه نقاط موجود در روش 
اجزای محدود استاندارد می‌باشد. ۳ مجموعة نقاطی 
هبشتن. که کر که از المان‌های آنها عبرر. کرفه‌است: و 
دراصطلاح آنها را بریده‌است. این نقاط در شکل (۲) با 


دایره مشخص شده‌است. ۱۳ مجموعة نقاطی است که 
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نوک ترک در المان‌های آنها قرار گرفته‌است این نقاط نیز 
در شکل (۲) با علامت مربع مشخص گردیده‌است [3]. 

همچنین در رابطهٌ (۱) ز توابع شکل استاندارد در 
هویساید می‌باشد. تابع هویساید در نقاط بالای ترک 
مقدار ۱+ و در نقاط پایین ترک مقدار 2-۱ را به خود 
می‌گیرد. تابع »9 نیز از حل تحلیلی الاستیک نوک ترک 
نشثت گرفته‌است و به‌صورت زیر بیان می‌ گردد: 


۳( 0 ۳ < [,ظ] 


نی ۹ 


شکل ۲ نقاط تحت تأثیر توابع عادی غنی‌سازی و 


توابع نوک ترک 


به‌دلیل وجود ناپیوستگی توابع هویساید که در 
ماتریس شکل ایجاد می‌گردند. انتگرال گیری به‌روش‌های 
عادی, نظیر استفاده از روش گاوس کارساز نیست و دقت 
لازم در انتگرال‌گیری وجود ندارد. در این حصوص 
روش‌های متعددی پيشنهاد شده‌است که یکی از آنها 
روش پیشنهادی دالبو می‌باشد. در این روش المانی که از 
آن ترک عبور کرده‌است به زیر المان‌های مثلثی تقسیم 
می‌شود و انتگرال گیری براساس این زیرالمان‌ها محاسبه 
می‌گردد. شکل (۳) نمونه‌ای از تقسیم‌بندی این 
زیرالمان‌ها را به تصویر می‌کشد [۵]. 
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شکل ۳ زیرالمان مثلئی منطقهُ ترک 


لته ور المان‌هایی که متطقة السان‌های تن راید 
المان‌های غنی‌شده مرتبط می‌سازند و به آنها المان‌های 
مخلوط (5اع«ع۱ظ عصن«316) گفته می‌شود. مشکلاتی 
ور قسگرانی تحرانی‌ها ماهاه یقرت که این مش کادارت 
را می‌توان با اصلاح توابع غنی‌سازی کاهش داد. به 
این‌صورت که توابع غنی‌سازی در المان‌های استاندارد 
صفر می‌شوند. در المان‌هایی که هم گره‌های آنها 
غت‌سازی, مس شو تلم بدون تقییر .هی مانقت و در المان‌های 


روش احزای محدود تطابقی 
روش اجزای محدود تطابقی به‌عنوان یک تکنیک موفق 
در آنالیز مسائل ترک به کار می‌رود. اساس کار این روش 
اصلاح مش و معرفی مش‌بندی مناسب براساس شرایط 
و نوع مسئله می‌باشد. اصولا در مسائلی که ترک وجود 
دارد به‌دلیل تمرکز بالای تنش در اطراف نوک ترک خطا 
در آن محدوده زیاد می‌باشد و لذا پاسخ مسئله به‌خصوص 
فاکتور شدت تنش خطای نسبتا بالایی خواهد داشت و با 
مقدار تحلیلی خود فاصله دارد. بر همین اساس روش 
اجزای محدود تطابقی با اصلاح مش و یافتن نقاط 
حساس مسئله که نیاز به ریزسازی بیشتر دارند. مقدار 
خطای مسئله را کاهش می‌دهد و درنتیجه جواب‌ها به 
مقدار تحلیلی آن نزدیک‌تر می‌گردد. معمولا از دو نوع 
رویکرد برای اصلاح مش در روش اجزای محدود تطابقی 
استفاده می‌شود: ۱- روش‌های مبتنی بر باقی‌مانده‌ها 
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(0005ع۱ 3960 اقبالهع1) و ۲- روش‌های مبتنی بر 
بهبود جواب‌ها (عمطع۱ 2560( «عه۲معع. در 
روش‌های مبتنی بر باقی‌مانده‌ها که درابتدا توسط بابوشکا 
و رینبولت ارائه گردید[6], نتایم حل به‌روش اجزای 
محدود در معادلات حاکم بر سیستم قرار داده می‌شود و 
خطا به‌صورت نرم انرژی و با محاسبةٌ باقی‌مانده‌های 
محلی این معادلات برروی مجموعه‌ای از المان‌ها 
(۵161) محاسبه می‌شود. از سوی دیگر در روش‌های 
مبتنی بر بهبود جواب که نخستین بار توسط زینکوویچ 
ی بان گردنی اوه کدف یرای «اطلاعات ستاو 
می‌شود تا مقادیر دقیق‌تری برای متغیرها به‌دست آید[7]. 
در این مقاله از این نوع رویکرد برای اصلاح مش استفاده 
شده‌است. 

در این روش برای ارزیابی خطاء بایستی از میزان 
اختلافی که بین جواب‌های بهبودیافته و جواب‌های قبلی 
وجود دارد استفاده می‌گردد. برای این کار ابتدا مسئله را 
با همان مش‌بندی اولیه آنالیز می‌کنيم. برای بهبود 
جواب‌ها کافیست بااستفاده از روابط توابع شکلی جواب 
بهبود یافته را محاسبه نماییم. این رابطه به‌صورت زیر بیان 
راهم ۸ 
(۳( 6( 6 

که "0 تنش بهبودیافته, تابع ااشکل و 6 تا 
تنش در گره‌ها می‌باشد. 

از آن‌جایی‌که جواب‌های بهبود یافته دارای تقریب 
مناسبتری نسبت به جواب‌های واقعی مسئله می‌باشد لذا 
می توان از اختلاف دو جواب به‌عنوان خطای تقریب 
مسئله استفاده نمود. 

(6) ۵ ده 

6جواب روش اجزای محدود در گره‌ها می‌باشد. 
تعریف فوق مقدار خطای برآوردشده را در تک‌تک نقاط 
دامنه می‌دهد. ولی برای ریزسازی شبکه معیار مناسبی 
نمی‌باشده چرا که ممکن است در یک نقطةٌ خاص مثل 
نوک ترک تنش به‌سمت بی‌نهایت برود که خطا در این 


نقطه بسیار بالا خواهد بو ولی درمجموع خطای کل 
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جواب قابل قبول باشد. ازاین‌رو به‌جای خطای 
نقطه‌به نقطه نرم خطا تعریف می‌شود که به‌صورت انتگرال 
عددی از تابع خطا روی دامنه است. یکی از معروف‌ترین 
نرم‌های خطاء نرم م1 می‌باشد که به‌صورت زیر تعریف 
قی گر هدز 
1/2 
(0( 0 ۱ -م 
9 


به این ترتیب نرم تابع خطای تنش برابر خواهد بود با: 


0( ظ_ - 6 7 ال ( -4- | -]مع| 


پس از تخمین خطای مسئله می‌توان یک شبکه بهینه 
برمبنای خطای برآوردشده ایجاد نمود. بدین ترتیب که 
در نقاط با خطای بالاتر المان‌های ریزتر و در مناطق با 
خطای پایین از المان‌های درشت‌تر استفاده نمود. با 
تعریف یک خطای هدف مطلوب ,( |,6) می‌توان 
تراکم المان‌های متصل به گره‌های مختلف را تعیین نمود 


و شبکة جدید را تعیین نمود. 
,6) 
- ۰( ظ) 


2۳ 
0 #ر حلا! 


0 
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برآورد فاکتور شدت تنش با ترکیب روش اجزای 
محد‌ود بسط یافته و تطابقی 


پارامتر فاکتور شدت برای بیان میزان تمرکز تنش در 
اطراف نوک ترک به‌کار می‌رود و نقش تعیین کننده‌ای در 
مسیر رشد ترک دارد. مفهوم فاکتور شدت تنش برای 
اتتشتیار تسط آروی تال ۱۵۷ ترا سشضی قدار 
تکینگی استفاده شد[9]. برای محاسبة فاکتور شدت تنش 
به‌صورت تحلیلی بایستی از روابط الاستیک موجود 
تاک یره که قفا هکره و مایا اس سا 
شده‌است و در ساير مسائل عموما از روش‌های عددی 
برای برآورد ضریب شدت تنش استفاده می‌شود. از 
روش‌های متداول عددی برای برآورد این ضریب می‌توان 
ب‌روش همبستگی تغییرمکان‌ها 


حعصص‌هاموزر۲] 


سال سی و یکم. شمارهُ سوم ۱۳۹۷ 


۵ «0۳612110)). روش گسترش ترک مجازی 
(0مطامجور صمتقصمان عفن لفتتاتز 6۷ روش انشگر ال 
بسته ترک اصلاح‌شده (متاومله علمفت .۰ ۷۲0۵01160[ 
موعامن؛_ روش [-انتگرال. روش انتگرال متقابل 
(0مصاعصط 1تععاصا مناعه۲ع1ع۲) اشاره داشت[2]. در این 
مقاله از روش انتگرال متقابل برای برآورد ضریب شدت 
تلقن انگفاده قنده‌امتت: 

روش انتگرال متقابل یکی از دقیق‌ترین روش‌ها 
تداع مان اک نیرسن اب رف 
به‌نوعی روش 1 انتگرال را ارتقا می‌دهد و بااستفاده از یک 
فان کمک اتف بارش تفت و مادک 
محاسبه نماید. 

انتگرال اندرکنش متقابل ۷1 برای اجسام الاستیک 
توتبنط باقو ی قزر سال: ۲۰۶۵ آرانه کردید[10] فش که 
این روش نسبت به سایر روش‌ها دارد این است که 
همزمان و با یک انتگرال‌گیری پارامتر فاکتور شدت تتش 
در هر دو مود محاسبه می‌گردد. در این روش بایستی یک 
فضای کمکی متأثر از تتش و کرنش در مسئله تعریف 
نمود. این فضای کمکی به‌نحوی در مسئله تعیین می‌گردد 
که معادلات تعادل و شرایط مرزی بدون نیرو برروی 
سطح ترک در منطقة ۸ را ارضا نماید. ۵۳ مساحت 
معادلی است که در آن انتگرال گرفته می‌شود. در این 
روش انتگرال 7 از ترکیب انتگرال مسئلة اصلی و مبثلة 


کشکی ار اشطرال انت‌رکتن وتا دم ابا 
0/5 + 1 ِِ 1 


حملات ۱ مربوط به ود ضعیت واقعی و کمک 
است و ]1۷ نیز انتگرال متقابل می‌باشد که مقدار آن برابر 


.10 ۲ 
۷ و ِ ]گم ۷ 


1 لا 1 1۳ [ 


(٩) 
برآع ار تباطاشکرال انز کیش ی ضریبا شلات تشن‎ 


از تعریف ضریب زیر استفاده می‌شود[10]. 


سال سی و یکم» شمارة سومء ۱۳۹۷ 


۶۳ 


2 
۰۱ دوک 0 - - ۷۲ 


برای مد 1 1< "5۳ و 0 سک فراومن کیزق و تزا 
مدلگل 0 15۳ و 1< 157 قرار می‌گیرد؛ بنابراین: 
:1 


1 - ۲ ۱۱) 
2 


ای "۳ برای حالت کرندٌ از و 
به‌جای ظ بر کرنش مسطح ی 
برای تنش مسطح نیز از ظ استفاده می‌شود. 

در اين برآورد مدل‌سازی ترک با روش اجزای 
محدود بسطیافته صورت می‌گیرد و هندسهٌ ترک از میان 
المان‌ها عبور می‌کند. پس از برآورد تنش‌ها و پیش از 
محاسبهٌ ضریب شدت تنش. خطای دامنه به‌روش اجزای 
محدود تطابقی تخمین زده می‌شود و درصورتی که خحطای 
دامنه بیش از محد‌وده قابل‌قبول باشد. شبکه با روش 
ارائه‌شده در بند (۲) محجدد تولید می‌شود و آنالیز شبکه 
می‌تواند از بین المان‌ها عبور نماید. زمانی‌که خطای 
تخمین زده‌شده در حد قابل قبول رسیك می‌توان بااستفاده 
از روش انتگرال متقابل ضریب شدت تنش را محاسبه 
نمود. با این ترکیب دیگر نیازی به دنبال کردن مسیر ترک 
تنش برمبنای تنش‌های با دقت قابل قبول انجام می‌گیرد. 

در تمامی مراحل آنالیزه شبکٌ اجزای محدود بدون 
توجه به مسیر ترک بدون تغییر بافی می‌ماند تا بتوان از 
مزایای روش اجزای محدود بسطیافته استفاده نمود. تا 
زمانی که خطای شبکه از حد قابل قبول بیشتر شود و تولید 
مجدد شبکه انجام گیرد. به‌اين ترتیب ممکن است در ده 
مره وف که فا بت مره ولد تاه 
گردد. شبکهة تولیدی جدید بدون نیاز به تطابق با مسیر 


ترک تولید می‌شود. 


مدل‌سازی عددی 


نشريهةُ مهندسی عمران فردوسی 


۶۴ 


معرفی‌شده دو مثال عددی مورد بررسی قرار گرفته‌اند. 
این دو مثال از مسائل کلاسیک مکانیک شکست هستند 
و ضرایب شدت تنش آن به‌صورت تحلیلی نیز موجود 
می‌باشد. برای آنالیز اجزای محدود از المان چهارگوش 
چهارگرهی استفاده شده‌است و انتگرال گیری اصلاح شده 
با چهار نقطةٌ گاوس صورت گرفته‌است. فرایند اصلاح 
ی لاسام .ار اواشله.کی شل. ۲ وت 
گرفته‌است. همچنین خطای هدف در هر دو مسئله برابر 
۵ درصد درنظر گرفته شده‌است. مثال نخست. مسئلهً 
ترک لبه‌ای تحت کشش و مثال دوم ترک میانی 1۵ درجه 
تحت کشش در نظر گرفته شده‌است. در مثال اول تنها 
مود اول فعال شده‌است اما مثال دوم مسئلة ترکیب مودها 
می‌باشد و کاوایی این تکنیک را در زمان فعال شدن هر 
دو مود شکست نشان می‌دهد. 


ورق مستطیلی با ترک لبه‌ای 
در این مثال یک ورق مستطیلی به‌طول 257 < 9 و ارتفاع 
000 2 20 مدل‌سازی شده‌است که دارای یک ترک 
لبه‌ای به‌طول اولیة 80 -ه می‌باشد (شکل 4). نمونة 
موردنظر از دو طرف با تنش یکنواخت که 
کشیده می‌شود. برای رشد ترک نیز در هر مرحله یک گام 
رشد ترک به‌طول 2 درنظر گرفته شده‌است. 
پارامترهای مکانیکی نظیر مدول الاستیسیته و ضریب 
پواسون» به‌ترتیب برابر ر 2-00 و ۷-03 
11 


۸+۸ ۸ 4 


تلاوت 


اعچعه 


07 


شکل ۶ ورق مستطیلی با ترک لبه‌ای 


نشرية مهندسی عمران فردوسی 


بهبود روش اجزای محدود بسط‌یافته د رآنالیز. 


برای نشان دادن کارایی تکنیک ارائه‌شده مسئله یک مرتبه 
تنها با روش اجزای محدود بسطیافته و بدون ترکیب با 
روش تطابقی انجام شده‌است و بار دوم در ترکیب دو 
روش صورت پذیرفته‌است. شبکٌ اولیه‌ای که برای آنالیز 
اجزای محدود بسط یافته ستفاده شده‌است» در شکل (6 
نشان داده شده‌است. در این شکل مشخص ات که ترک 
از بین المان‌ها عبور کرده‌است. 


0 ۱ 
بو 


3 
0 
3 
ِ 
3 
0 
3 
3 


شکل ٩‏ شبکه اجزای محدود بسط یافتةً اولیه 


معیاری که برای رشد ترک در هر دو فاز درنظر 
گرفته شده‌است براساس معیار حداکثر تنش محیطی 
شدت تنش دو مود شکست به‌صورت زیر محاسبه 


یه ۳ 


2 

)۱۲ 8+ 1 ب | ]هه 
4۱ »41 21 

تخمین خطا در هر مرحله از رشد ترک صورت 
گرفته‌است ولی اصلاح شبکه‌ای انجام نشده‌است و 
درنتیجه با رشد ترک میزان خطا افزایش يافته و از مقدار 
خطای هدف ۱۵ درصد نیز فراتر رفته‌است که نمودار 
تغییرات آن در شکل (0) آورده شده‌است. علت این امر 
حرکت نقاط حساس مسئله (منطقةٌ اطراف ترک) با رشد 
ترک می‌باشد. چون منطقهٌ ریزشده قبلی دیگر با نقاط 
حساس جدید مسئله منطبق نیستند خطای حل شروع به 
افزایش می‌نماید که با اصلاح مجدد شبکه این مسئله 


سال سی و یکم. شمارهٌ سوم ۱۳۹۷ 


برطرف می‌شود. 


14 16 18 


12 


طول ترک (سانتیمتر) 


10 


درصد خطا 


شکل 1 تغییرات خطا با رشد ترک بدون اصلاح شبکه 


سال سی و یکم شمارُ سوم ۱۳۹۷ 


۶۵ 


بااستفاده از روش انتگرال متقابل محاسبه شده و در 
جدول (۱) آورده شده‌است. همچنین این مقدار با مقدار 


تحلیلی این مود که از رابطة زیر به‌دست می‌آید[11] 


مقایسه شده‌است. 


16, 2- 2 


و 


(۳ 


بدیهی است که در این مسئله مود دوم فعال نشده‌است و 


قیرفته شلات ی اضق اسخه 


درصد خطا | 14 تحلیلی 


1۲ 
۳ 
۱۳/3۹ 
۳2 
۲۱۸۳۷ 


۲۳,۲ 


۸۱۳ 
۱۱۹/۷۸ 
۱,۷ 
۲۳۹۸۹ 


۳۳,۹ 


2۳,۳۹ 


4 عددی 


ال عددی مختصات نوک ترک 
۳2 (۶۹,۹۰و )٩,۹۹‏ 
۰ - (۶۹,4۳و۱۱,۹) 
۵ - (۵۰,۰۶و۱۳,۹۸) 
۱,1۹4 (۶۹,۹۱و۱۵,۹۸) 
۱ (۶۹,۲۱و۱۷,۸۵) 
۳,۵- (1۸,۸۸ و ۱۹,۸۲) 


گام ۶ اصلاح یافته گام ۴ اصلاح یافته تاء گام ! الی ۴ 


شکل ۷ شبکة اصلاح‌شده در طول رشد ترک 


نشریهُ مهندسی عمران فردوسی 


۶۶ 


حال در مرحلهٌ دوم ترکیب روش‌های اجزای 
محدود بسطیافته و تطابقی انجام گرفت و ملاحظه گردید 
که در دو مرحله از رشد ترک. خطا از حد مجاز ۱۵ درصد 
فراتر رفت و نیاز به اصلاح شبکه وجود دارد. شبکةٌ 
اصلاح‌شده در این دو مرحله در شکل (۷) نشان داده 
شده‌است که تراکم المان‌ها در نوک ترک هرمرحله, نشان 
می‌دهد روش تطابقی به‌خوبی نقاط حساس مسئله را 
یافته و شبکة مناسبی تولید کرده‌است. 

میزان تغیبرات خطا در مراحل مختلف رشد ترک 
در شکل (۸) آورده شده‌است که ملاحظه می‌شود در 
مراحلی که خطا بالای خطای هدف است. اصلاح شبکه 
منجر به افت خطا شده‌است. حال با روش ترکیبی مجددا 
برآورد ضریب شدت تنش انجام شده‌است که نتایج و 
مقایسٌ آن با مقادیر تحلیلی در جدول (۲) نشان داده 


مه اببتتتا, 


خطای هدفه 


درصد خطا 


طول ترک (سانتیمتر) 


شک ۸ تغییرات خطا با رشد ترک با روش ترکیبی 


بهبود روش اجزای محدود بسط یافته د رآنالیز. 


با مقايسة مقادیر جداول (۱ و ۲) واضح است که 
استفاده از روش ترکیبی به‌مقدار قابل‌توجهی خطای 
آنالیز مجدد روی تمامی شبکه‌ها انجام گرفته‌است و 
نیازی به فرایند انتقال اطلاعات نمی‌باشد. 


ورق مربعی با ترک مورب 

مدل‌سازی دوم عددی مورد بررسی ترک مورب ۶۵درجه 
در داخل یک ورق مربعی می‌باشد که تحت کشش قرار 
قایف انم مسظله ه ای اط اشت دایه که یک وا 
مود ترکیبی است که هر دو مود در آن فعال می‌شوند و 
کاواس فکن که اراه‌شته در این کرته سانشان 
می‌دهد. در اين مثال عرض صفحه 0-506 و ارتفاع 
50 < 20 درنظر گرفته می‌شود. طول اولية ترک نیز 
104/20 < ۵ می‌باشد. نمونهُ مورد نظر از دو طرف با 
تتش یکنواخت برابر ,1۳8و کشيده می‌شود. برای 
رشد ترک نیز در هر مرحله یک گام با طول 2 در نظر 
گرفته شده‌است. ساير پارامترهای مکانیکی نظیر مدول 
الاستیسیته و ضریب پواسون به‌ترتیب برابر 
م۱000 - 8 و ۷۲0.3 می‌باشند (شکل 4). 


جدول ۲ مقادیر ضرایب شدت تنش با ترکیب روش‌ها 


مرحله درصد خطا | أع تحلیلی 
اول ۷:۷ ۸٩۱‏ 
دوم ۱۰,۵۱ ۱۱۹/۸ 
سوم ۱۳۹ ۱۷ 
چهارم ۱2۷۶ ۲۱۳۸۹۸۹ 
چهارم اصلاح‌شده ۲« ۲۲۵۹۸۵ 
پنجم ۱۳۹۵ ۳۳,۲۳۸ 
ششم ۱۳5۵ 1۸,۰ 
ششم اصلاح‌شده ۹,۲۰ 1۸,۰ 


نشرية مهندسی عمران فردوسی 


6 عددی | ]61 عددی | مختصات نوک ترک 
۱۷,9۶ ۳ (4۰,٩۶و۸۹,٩)‏ 
۹/۳۹ ۸ - (۳, ٩۶و‏ ۱۱,۹) 
۱۳/۳۸۹ ۲ (۵۰,۰۶و۱۳,۹۸) 

۱۸۳۰۸ ۱2۶ ِ 
۳ ۱۳ ۰:۱۷ (۶, ۵۰و ۱۵,۹۸) 
۱۳۰/۳۶ ۱,۱۸ (۶۹,۸۵و۱۷,۹۸) 

۹ ۱۳۸۰۵ 
)۱۹:۹۸۵ ۶۹,۸۸۱ ۳ 329 


سال سی و یکم. شمارهٌ سوم ۱۳۹۷ 


۲۷۸۸۱ 


شکل ٩‏ ورق مربعی با ترک مورب 


مطابق مثال قبل این مسئله نیز به دو صورت آنالیز 
می‌شود؛ بار نخست با اعمال روش اجزای محدود 
بسط یافتةُ تنها و مرحلة دوم با اعمال روش ترکیبی. از 
آن‌جایی که مقدار تحلیلی ضریب شدت تنش فقط در 
مرحله اول و قبل از رشد ترک موجود می‌باشد لذا برای 
بیان هم‌گرایی ضریب شدت تنش. گام اول یک بار با مش 
اولیه و بار دوم با مش اصلاح‌شده حل می‌شود که این دو 
شبکه در شکل‌های (۱۰ و ۱۱) نشان داده شده‌اند. 


شکل ۱۱ شبکة اجزای محدود بسطیافته اصلاح شده 


مقادیر ضرایب شدت تنش بااستفاده از هر دو شبکه 
محاسبه شده‌است و با مقدار تحلیلی آن که از رابطة (۱۶) 


نشریة مهندسی عمران فردوسی 


۶۷ 


به‌دست می‌آید [11] مقایسه شده‌است و در جداول (۳ و 
۶ آورده شده‌است. باتوجه به وجود دو نوک ترک در 
مسئله این برآورد برای هر دو نوک ترک ابتدایی و انتهایی 
انجام شده‌است. 

13 معلاهک0 > ررک < رک 


تغیاول ۳ مقادیر ضرافب شلات ففی با شک ارلیه 


نوک ابتدای ترک | نوک انتهای ترک 
| لک | ۲ ] لک | ۲ ۲ لک 


۳۶۸ ۳۰/۱ ۳۵۶ ۳۶۲ ۳۳۵ ۳۳۵ 


جدول ۶ مقادیر ضرایب شدت تنش با شبکة اصلاح‌شده 


نوک ابتدای ترک | نوک انتهای ترک 
1 1 آع 1 1 1 


۲۵۰ | ۲۸ ۲۳۸ | ۲,۵ | ۲۳۵ ۳۵ 


خطای هدفه 
خطای مساله و 20 
1 15 وه 
1 
۱ 
3 
10 3 
۳۷ 
۱ رخ 
40 30 20 10 0 
طول ترک (سانتیمتر) 


شکل ۱۲ تغییرات خطا با رشد ترک بدون اصلاح شبکه 


نتایج دو جدول حاکی از افزایش دقت پس از 
اعمال اصلاح شبکه می‌باشد. پس از محاسبة ضرایب 
شدت تنش فرایند رشد ترک انجام گردید که مشاهده شد 
ترک از حالت مورب شروع به چرخش به‌صورت مستقیم 
کرد. فرایند رشد ترک با روش اجزای محدود بسط یافته 
تنها نشان‌دهندة افزایش خطا در مراحل مختلف رشد ترک 
می‌باشد که در شکل (۱۲) نشان داده شده‌است. 


سال سی و یکم. شمارهُ سوم ۱۳۹۷ 


۶۸ بهبود روش اجزای محدود بسطیافته د رآنالیز. 


درحالی که با ترکیب روش با اجزای محدود تطابقی در در روش اجزای محدود بسطیافته کنترلی برروی میزان 
دو مرحله (گام‌های ۶ و 7 رشد ترک) مقدار خطای مسئله خطای مسئله وحود ندارد. و در روش اجزای محدود 
در محدودهٌ مجاز نگه داشته شده‌است که در شکل (۱۳) تطابقی با هر گام از رشد ترک شبکه متناسب با مسیر 
نمایش داده شده‌است. ترک باید دوباره از ابتدا تشکیل شود. با ترکیب این دو 
روش می‌توان به‌شدت این مشکلات را کاهش داد. به این 

خطای ماه وید انجام شود و ترک از داخل المان‌ها عبور کند و دیگر با 

خطای حل بیش از حد قابل قبول گردده فرایند اصلاح 
10 مش صورت پذیرد. این تولید مجدد شبکه کاملا مستقل 
1 24 ِ ِ ِ 0 مرحله از رشد ترک. مسئله با یک شبکة یکسان حل شود. 


درصد خطا 
۵ 


طول ترک (سانتیمتر) نتایج مدل‌سازی‌های عددی نشان دادند که درصورت 
ترکیب این روش‌ها. خطای حل به‌میزان قابل توجهی 


کاهش می‌یابد و برآورد ضرایب شدت تنش با دقت 


شکل ۱۳ تغییرات خطا با رشد ترک با روش ترکیه 
تیور بر ان 


ی بیشتری انجام خواهد شد. این تکنیک در ترک‌های 
نتیجه گیری ۱ ۲ ِ ۲ 

تک‌مود و هم در مودهای ترکیبی قابلیت برآورد دقیق را 
روش اجزای محدود بسطیافته و روش اجزای محدود 
تطابقی هر کدام قابلیت‌های بالایی برای حل مسائل آنالیز 


نشان داده‌است. 
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